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Abstract. We perform a systematic WKB expansion to all orders for a 
one-dimensional system with potential V(x) — Uq/ cos 2 (ax). We are able 
to sum the series to the exact energy spectrum. Then we show that at any 
finite order the error of the WKB approximation measured in the natural 
units of the mean energy level spacing does not go to zero when the quantum 
number goes to infinity. Therefore we make the general conclusion that 
the semiclassical approximations fail to predict the individual energy levels 
within a vanishing fraction of the mean energy level spacing. 
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In the last years many studies have been devoted to the transition from 
classical mechanics to quantum mechanics. These studies are motivated by 
the so-called quantum chaos (see Ozorio de Almeida 1990, Gutzwiller 1990, 
Casati and Chirikov 1995). An important aspect is the semiclassical quan- 
tization formula of the energy levels for integrable and quasi-integrable sys- 
tems, i.e. the torus quantization initiated by Einstein (1917) and completed 
by Maslov (1972, 1981). As is well known, the torus quantization is just 
the first term of a certain ^-expansion, the so-called WKB expansion, whose 
higher terms can be calculated with a recursion formula at least for one degree 
systems (Dunham 1932, Bender, Olaussen and Wang 1977, Voros 1983). 

Recently it has been observed by Prosen and Robnik (1993) and also 
Graffi, Manfredi and Salasnich (1994) that the leading-order semiclassical 
approximation fails to predict the individual energy levels within a vanishing 
fraction of the mean energy level spacing. This result has been shown to 
be true also for the leading (torus) semiclassical approximation by Salasnich 
and Robnik (1996). 

In this paper we analyze a simple one-dimensional system for which we 
are able to perform a systematic WKB expansion to all orders resulting in 
a convergent series whose sum is identical to the exact spectrum. For this 
system we show that any finite order WKB (semiclassical) approximation 
fails to predict the individual energy levels within a vanishing fraction of the 
mean energy level spacing. 

The Hamiltonian of the system is given by 

H=£- + V(x), (1) 

where 

Vtx) = ■ (2) 

cos (ax) 

Of course, the Hamiltonian is a constant of motion, whose value is equal 
to the total energy E. To perform the torus quantization it is necessary to 
introduce the action variable 

1 r \/9m i — 

/ = I pdx = ^(VE - Ju ) . (3) 

The Hamiltonian as a function of the action reads 
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and after the torus quantization 



I = (v + \)h, (5) 



where v — 0, 1, 2, . . ., the energy spectrum is given by 



EF = A[( V + \) + \b\\ (6) 

where A = a 2 h 2 /(2m) and B = ^8mUo/(ah). 
The Schrodinger equation of the system 

l-^^ 2 +V(x)Mx) = E^x), (7) 

can be solved analytically (as shown in Landau and Lifshitz 1973, Fliigge 
1971) and the exact energy spectrum is: 

El* = A[{v+\)+ l -Vl^W]\ (8) 

where v = 0, 1, 2, . . .. We see that the torus quantization does not give the 
correct energy spectrum, but it is well known that the torus quantization is 
just the first term of the WKB expansion. To calculate all the terms of the 
WKB expansion we observe that the wave function can always be written as 

i 

ijj(x) = exp (-o-(rr)) , (9) 

where the phase a(x) is a complex function that satisfies the differential 
equation 

a'\x) + (-)a"(x) = 2m(E - V(x)) . (10) 
The WKB expansion for the phase is given by 

fc=0 1 

Substituting (11) into (10) and comparing like powers of H gives the recursion 
relation (n > 0) 

a'l = 2m(E - V(x)), £ a' k a' n ^ k + a" n ^ = . (12) 

k=0 
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The quantization condition is obtained by requiring the single- valuedness 
of the wave function 



/°° Ti r 
da = ]T(-) fc f da k = 2nhv , (13) 



where v — 0, 1, 2, . . . is the quantum number. 

The zero order term, which gives the Bohr-Sommerfeld formula, is given 

by 

<fda = 2 J dx^2m(E - V(x)) = 2irh(^ - ^B) , (14) 

and the first odd term in the series gives the Maslov corrections (Maslov 
index is equal to 2) 



! j) j> dcr l = ( j) J ^Plcontour = • (15) 



The zero and first order terms give the equation (6), which is the torus 
quantization formula for the energy levels (Bohr-Sommerfeld-Maslov). Here 
we want to analyze the quantum corrections to this formula. We observe 
that all the other odd terms vanish when integrated along the closed contour 
because they are exact differentials (Bender, Olaussen and Wang 1977). So 
the quantization condition (13) can be written 



00 h r 1 

E(-) 2 " f d ^k = 2vr% + -) , (16) 
fc=0 1 J 1 



thus again a sum over even-numbered terms only. The next two non-zero 
terms are (Narimanov 1995, Bender, Olaussen and Wang 1977, Robnik and 
Salasnich 1996) 

(-) 9 da 2 = -—=-—- / dx = \ > ■ , (17) 



A*/, ^ 1 d 3 r V"\x) 1 d 4 r V' 2 (x)V"(x) 1 

V f da4 - (2^)3/^120^7^^-^) 288dE*J dX ^V{x) 1 ' 

(18) 
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A straightforward calculation of these terms gives (see the Appendix) 

2nh 



( )»^= (19) 



and 



h 4 /" , 27r/i 



Up to the fourth order in h ~ -B 1 the quantization condition reads 

EW=A\(i/+-) + -B + — — l 2 . (21) 

2 2 AB 16fi 3j 1 J 

The first two terms on the right side give the torus quantization formula, 
and the other two terms are quantum corrections. Higher-order quantum 
corrections quickly increase in complexity but in this specific case they can 
be calculated. We first verify by induction, following Bender, Olaussen and 
Wang (1977), that the solution to (12) has the general form 

a' n (x) = {a' y- 3n P n {cos (ax)) sin^ (ax) , (22) 

where f(n) — for n even and f(n) = 1 for n odd, and P n is a polynomial 
given by 

P„(cos (ax)) = ]T C n>l cos 2 '" 3 " (ax) , (23) 

1=0 

with g(n) = (3n — 2)/2 for n even and g(n) = (3n — 3)/2 for n odd. 

The integrals in (16) are performed by substituting z = tan (ax). In this 

way the 2/c-term reduces to 

~i J = <7> f C ^t dz (E-U -U 0Z ^ 

(24) 

We observe that 

(1 + z 2 ) 3k-l-l 



dz 



(E-U - Uoz 2 )^- 1 / 2 



,3,-i r(i) L ^ (i + z *r-^ 



( 1} ~ r(3k- \) dE**- 1 f dz (E-u - u z*y/i ' (25) 
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so the only non-zero term is for I = 

d 3k-i (l-j-^afc-i ^ 2 6fc ~ 1 T(3A; - 1/2) 2 d 3 ^ 1 3fc _ x 



dE^~ l J (E-U - U z 2 y/ 2 c/V2 r(6A; - 1) <9£ 3fc 

r(3fe-l/2) 2 1 

= ttf' 2 r(6fc-i) r ^^72 27r ' ^ 

where (5 — {E — U )/U . At this stage we obtain 

(5)»/<fo* = (_i)--. R -P!!!^ C2M _L_ 2;r . (27) 

Now we need to find the coefficient C2fc,o explicitly. By inserting (22) with 
(23) in the recursion relation (f2) we obtain 

n n— 1 

C kfi C n - kfi -(2mU a)C n - lfi = £ C kfi C n - kfi +2C nfi -(2mU a)C n - lfi = , 

k=0 k=l 

(28) 

from which we have 

^ i 

C kja = -[(2maC/ )C fc _i i0 - £ C jt oC k - jfi ], C 0fi = 1 . (29) 

From this equation one shows C^o = maUo. Further, it easy to show that all 
higher odd coefficients vanish, C*2fc+3,o — for k — 0, 1, 2, . . .. The solution of 
this equation for the remaining nonzero even coefficients is given by 

C 2kfi =(-l) k (2mU a) 2k 2' 2k (^j , (30) 

which can be verified by direct substitution in equation (29) resulting in an 
identity for half integer binomial coefficients. Then the integral (27) can be 
written 



fda 2k = (-l)^W- 1 (2m)^-*2-*gj^-i/2 = -ig) . 

(31) 
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In conclusion, the WKB quantization to all orders (16) is 



i i 00 /IN i 
^ = A[( V+ -) + -^^f. (32) 

Because ET=o(l) Bl ~ 2k = Vl + B 2 we have E e u x = E^°\ i.e. the WKB 
series converges to the exact result (8). 

Now we can calculate the error in units of the mean level spacing AE U = 
El x +1 — E ex between the exact level E ex and its WKB approximation E^ 
to Nth. order: 

" — = - V \A-±—, forz/^oo. (33) 

The limit clearly shows that even for arbitrarily small but finite h (1 << B < 
oo), the relative error for any finite WKB approximation becomes constant 
on increasing u, and scales as 

E ex _ E {N) X / I N 1 

^ae^~2{nIi)b^' B ^™' < 34 » 

Note that the limit B — > oo is equivalent to the limit H — > 0. 

For our present system we can conclude that to any finite order semiclassical 
approximation the error measured in units of the mean level spacing remains 
constant even if the quantum number increases indefinitely, contrary to the 
naive expectation. This confirms the general statements made by Prosen 
and Robnik (1993). We have thus provided a clear demonstration that the 
semiclassical methods cannot predict the individual energy levels (and also 
their wavefunctions) within a vanishing fraction of the mean energy level 
spacing. Therefore we cannot expect the semiclassics to correctly describe 
the fine structure of energy spectra manifested in the short range statistics 
like the energy level repulsion, which was predicted to be a purely quantum 
effect (Robnik 1986), later reconfirmed by Berry (1991). On the other hand 
Prosen and Robnik (1993) have shown that the long range statistics of the 
energy spectra are very well captured even by the lowest order semiclassical 
approximation. This is of course compatible with the very important semi- 
classical theory of delta statistics A(L) (spectral rigidity) by Berry (1985), 
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employing the Gutzwiller periodic orbit theory (1990), where agreement with 
predictions of random matrix theories and with the experimental and numer- 
ical data has been obtained at large L. Also, Berry and Tabor (1977) have 
used torus quantization of integrable systems (with many degrres of free- 
dom), predicting the Poissonian (exponential) energy level distribution. Our 
results show that their result cannot be rigorous, especially as we know some 
counterexamples of integrable systems with non-Poissonian statistics (Bleher 
et al 1993), and also know that their approximation does not take into ac- 
count the nonperturbative tunneling effects, but it is nevertheless a heuristic 
argument explaining why typically we do observe Poissonian statistics in clas- 
sically integrable systems. By typically we mean that the set of exceptions 
has a small or even vanishing measure. 

The conclusion of this paper is that the semiclassical methods are just not 
good enough (at any order) to describe the fine structure of energy spectra 
and wavef unctions. Our approach leading to the above conclusion rests upon 
a systematic WKB expansion for the potential V(x) — Uq/ cos 2 (ax) using 
the technique of Bender, Olaussen and Wang (1977). We are able to calculate 
all orders, the series is convergent and can be summed precisely to the exact 
result. 
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Appendix 



In this appendix we show how to obtain the formulas (19) and (20). In all 
integrals of this section the limits of integration are between the two turning 
points. After substitution z = tan (ax), we have 

V'\x) AaU 2 ryfi z\z 2 + 1) 



f dx v _ m = r 



dz- 



where (3 — (E — U )/U . In conclusion we have 



Ui %f + 4P)1), (35) 



: i»l fc = -s« = -¥' (36) 



with B = y/8mU /(ah). 

To obtain the formula (21) we proceed in the same way. 

fdx V " 2 ^ 4a 3 £/ 2 ryfi ^ {9z 4 + 6z 2 + l){z 2 + 1) 

' x ^e- v(x) VUo~ z VP-? 2 

= -^(45/3 3 + 90/3 2 + 56/3 + 16)^, (37) 

From which we obtain 

m J iix ^wr ■■ (38) 

For the last integral we have 

V' 2 (x)V"(x) 8a 3 U 2 ry/P J z 2 (3z 2 + l)(z 2 + I) 2 

dz- 



J Jf, - VOA ^/C^ J- 



z 2 



SaU jl(105p 4 + 280/3 3 + 240/3 2 + 60/3) , (39) 



'U y ' 128 
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from which we obtain 
d 4 



In conclusion we have 



<7>' 



doA 



V ,2 {x)V"{x) _ 3157m 2 



dx 



dE'J p - V(x) 



h 4 



a 



TT a , 1 135 



J J 2 



(2m) 3 / 2 
64(2m) 3 /2[/2 



120 2 

2nh 



165 3 
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